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The emergence of advanced lightweight materials is resulting in a new generation of 
lighter, flexible, more-efficient airframes that are enabling concepts for active aeroelastic 
wing-shape control to achieve greater flight efficiency and increased safety margins. These 
elastically shaped aircraft concepts require non-traditional methods for large-scale multi- 
objective flight control that simultaneously seek to gain aerodynamic efficiency in terms of 
drag reduction while performing traditional command-tracking tasks as part of a complete 
guidance and navigation solution. This paper presents results from a preliminary study of a 
notional multi-objective control law for an aeroelastic flexible-wing aircraft controlled 
through distributed continuous leading and trailing edge control surface actuators. This 
preliminary study develops and analyzes a multi-objective control law derived from optimal 
linear quadratic methods on a longitudinal vehicle dynamics model with coupled aeroelastic 
dynamics. The controller tracks commanded attack-angle while minimizing drag and 
controlling wing twist and bend. This paper presents an overview of the elastic aircraft 
concept, outlines the coupled vehicle model, presents the preliminary control law formulation 
and implementation, presents results from simulation, provides analysis, and concludes by 
identifying possible future areas for research. 


I. Introduction 

C ontinuing advances in modem engineered materials and the need for greater efficiency of aircraft in flight are 
driving new aircraft design concepts that incorporate lightweight flexible multi-functional structures. These new 
structural concepts provide less rigidity while maintaining the same load-carrying capacity as previous generation 
designs. The Boeing 787 Dreamliner aircraft’s highly flexible wing structure is an example of this emerging trend. 
These advances are also driving research into active aeroelastic shape control of these flexible aircraft structures during 
flight. For instance, elastically shaping an aircraft by actively controlling wing wash-out twist and wing bending 
deflection in order to change the local angle of attack has the potential to yield lower fuel burn by drag reduction 
during cruise. Structural flexibility can further be leveraged to realize a revolutionary, optimal wing shape design that 
can accommodate a significant curvature for drag reduction benefits. 

A conceptual study was recently conducted outlining the concept for Adaptive Aeroelastic Shape Control 
(AASC) 1 . This study outlines plans for multi-disciplinary design, analysis and optimization to examine the potential 
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benefits of AASC on a flexible aircraft in flight over a conventional aircraft with traditional flight control system 
design. 

Four major technical areas of AASC research were outlined in this study: 

1. Vehicle concept design and optimization 

2. Aeroelastic flight dynamic modeling 

3. Elastically wing shaping actuation design 

4. Flight control design and vehicle simulation 

Significant research activity has occurred on the AASC in several of these areas, including a distributed actuation 
system for effecting elastic shape control 2 , investigation of non-linear coupling between aeroelastic modes and the 
flight vehicle dynamics 3 , development of a coupled high-fidelity numerical model utilizing vortex-lattice and finite- 
element methods 4 , offline shape optimization studies for drag minimization 5 ’ 6 , and development of a lower-order 
aeroelastic vehicle model that is appropriate for control law development 1 . 

This study builds on results from this previous research and focuses on flight control design in technical area 4 
above. Research under this area seeks to develop multi-objective flight control to simultaneously gain aerodynamic 
efficiency while maintaining traditional pilot command-tracking tasks for guidance and navigation. Research in this 
area also seeks to develop guidance laws to achieve low drag objectives, with specific focus on cruise segment of 
flight, applying modern multi-objective optimal control techniques for control of elastically-shaped aircraft. 

This paper presents a preliminary study of a notional multi-objective control law for an elastically-shaped flexible- 
wing aircraft with distributed leading and trailing edge control surface actuators. The purpose of this initial study is 
to determine feasibility of a centralized control framework that simultaneously minimizes drag, stabilizes structural 
modes, and simultaneously performs longitudinal attack-angle command-tracking tasks. This initial study continues 
the previous research by implementing the proposed control concept 7 , utilizing distributed control surface actuation 
for control of the coupled aeroelastic and flight dynamics model. This study develops a multi-objective control 
structure in a standard unconstrained linear-quadratic optimization framework and implements this controller on the 
recently developed model. The goals of this initial study are to demonstrate feasibility of this approach using a 
preliminary control concept, provide an assessment of expected performance improvements, identify challenges, and 
outline directions for future research to enable large-scale multi-objective control of aeroelastically shaped vehicles. 

This paper is organized as follows. The concept for elastically-shaped aircraft is presented, followed by a coupled 
aeroelastic flight dynamics model on which this work is based. The optimal control problem is formalized, the solution 
is presented, and the simulation implementation is presented. The closed-loop response of the system is presented 
with a focus on estimating quantifiable performance benefits. The results will be discussed, limitations identified, and 
a directions for future work are presented. 

II. Elastically Shaped Aircraft Model 

The model 4 utilized in this study describes the longitudinal dynamics of a generic transport aircraft at Mach 0.80 
and Mach 0.88 at an altitude of 35,000 ft. The longitudinal dynamics are described by two aircraft states: angle of 
attack a and pitch rate q. The aircraft model contains 80 structural states capturing the 20 dominant flexible modes, 
and provides 23 control inputs: one elevator, 11 flaps, and 11 slats; as shown in Figure 1. 
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Figure 1. Distributed leading-edge (slats) and trailing-edges (flaps) actuator configuration. 


The coupled equations of motion are given by 



( 1 ) 


where the states are given by x r — [ a q] T , x e — [q w qg q w qg] T , and q w ,q g £ are the generalized 
coordinate vectors for bending and torsion, respectively, along with their derivatives q w , q g £ M m . The elements of 
q w and q g are indexed such that q w . and q g . are the wing bending and torsional deflections of mode i = 1, . . ,,m, where 
m=20 is the total number of mode shapes. The input u r = [5 e ] is the elevator flap deflection angle, and u e = 

| by S s | 7 where Sf is the vector of 1 1 distributed flap inputs and S s are the 1 1 distributed slat input. This results in a 

model with 82 states and 23 inputs. This initial study assumes all aeroelastic modes and vehicle states are observable. 

The lift and drag characteristics of the aircraft are functions of the current vehicle states, flexible-mode states, and 
control inputs as described by the following equations 


AC d - K(2C l AC l + ACl ) + C Dy u + C d ^ 2 u 2 

C L = C L0 +AC L “ (2) 

AC l = C Lx x + C Ly u 


where K is the drag polar parameter, C D is the total drag coefficient, and C L is the total lift coefficient. 


III. Control Law Formulation 

The control system is designed using full-state feedback with the addition of tracking error integral feedback as 
developed in the literature 4 and summarized below. The controller tracks a reference angle of attack signal, a c . An 
integrator state x a integrates the attack angle error signal for feedback, given by 

x a A a (cr r rr c ) (3) 


where A a is the integrator gain, and can be combined with the system in (1) to form an augmented system for control 
development. The final augmented system is then of the form 


x — Ax + Bu + z 
y = Cx + Du 

where x £ is of size N=83 and u £ R M is of size M=23, given by 
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( 4 ) 


( 5 ) 
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The feedback control cost function J includes quadratic costs terms to address the multiple control objectives. 

( 6 ) 


J(x,u,t) = - (x T Qx + u T Ru + k d C D )dt 

Jo 

The weighting matrix Q = diag{Q r , Q e , Q a }, Q > 0 and its submatrices are diagonal matrices to achieve tracking 
of the vehicle states x r , control of the aeroelastic states x e , and the a command tracking integral, respectively. 
Similarly, R = diag{Rf,R s ,R e \,R > 0 and its submatrices penalize control actuation effort of the flaps, slats, and 
elevator, respectively. Total drag C D term with a scalar weight k d > 0. Substituting C D into J yields 


/O 


If 00 / 2 

, u, t) = — J [x T Qx + u T Ru + k d C Do + k d K(C Lo + C Lx x + C Lu u) + k d C Du i 

+ u T k D diag ( C D u) dt 


The control solution for problems of this form 7 is given by 

u(x,z ) = u x x + u z z + u c 

where 

u x = -R- 1 {B T P_+q d KCl u C Lx ) 
u z = -( R~ 1 B T S ) 


u c = —R 


- 1 (q d KC[ u C Lo +^q d CE u +B T A 0 ) 


The value of P is found from the solution to the following algebraic Riccati equation 

PA+A r P - PBR~ 1 B t P + Q = 0 

where 

R - R + q d KC[ u C Lu + q d C D ^ 2 
A=A- BR~\ d KCl u C Lx 
Q = Q + <7 d KC[ x (C Lx - C Lu R-'q d KCl u C Lx ) 


and 


5 = (PBR~ 1 B T -i 7 )- 1 ? 

A 0 = ( PBR~ 1 B t - i 7 )- 1 ( q d KC T Lx C Lo - ( PB + q d KC[ x C Lu )R ~ 1 (< l d KC T Lu C Lo +^q d C T Du 


(7) 


(8) 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 


The sufficient condition for a solution to exist requires checking the following, from which Q > 0 can be implied. 


k d KC L C Lx < R 


(13) 


IV. Implementation 

Forming the closed loop system of (4) under the feedback law (8), and noting that z — —A a va c where 
v = [0 ... 0 l] 7 , yields 

x = (A + Bu x )x — A a (Bu z + /)i7(a c ) 

u — u x x — A a u z va c + u c (14) 

y = x 

Without loss of generality the constant term Bu c has been dropped from the derivative equation in (14), noting 
that this term results in an offset x s of the state vector, which is given by 


x — x + x s 

x s = A-'BR- 1 (q d KC T Lu C Lo +^q d C T Du + B 7 T 0 ) 


( 15 ) 
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where x is the complete solution state vector. Also note that equation (15) can be solved using the original system 
\A rr A re 

matrix . . , which is invertible, and the effect of the constant u c term in the integrator can be solved separately. 

Aer ™ee. 

For evaluation of total drag, we define a total drag metric which integrates C D over the specified time interval. 

tf 


'avg 


if 

t/J 0 


= — | AC D dt 


(16) 


The state and aeroelastic weighting matrices were initially selected based on Bryson’s criteria, which sets the initial 
weighting elements based on the maximum expected value of the state and input elements. The Q and R diagonal 
matrices are given by 


Qe kqe 


Q r — kqr ' cLiag 
diag[l/qi max 1/qi 
Qa kqa 


R r k re 


(Y 2 

^ max 
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Qmax 


R e — diag j k r f * diag 
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v(u ei y 


e max- 

, k rs ’ diag 


1/ql 

I ~Wrt 


(17) 


maxex ; 


■p(u ei ) . 


The scalar parameter set k = {k d ,kq r ,kq e ,kq a ,k re ,k rsi k r f} allows for tuning of relative weights in the cost 
function as permissible, subject to the conditions in (13). 


V. Results 

Evaluation of the open loop models for both Mach 0.80 and Mach 0.88 shows that the two models are similar in 
structure and have unstable open-loop modes, but the both pairs (A, B) are controllable and the unstable modes can 
be stabilized through LQR feedback, as exemplified in Figure 2. 



Re Axis 


Figure 2. Open-loop and closed-loop poles. Shown for Mach 0.80 model under 
uniform feedback, the closed-loop system poles are all in the left-hand plane. 

A. Step-Response Study Results 

The model and controller were implemented in Matlab for evaluation. To evaluate the control effectiveness and 
drag reduction efficacy, the results were compared against a single-objective attack angle tracking controller, which 
was approximated by setting the state tracking error to a large value and setting the rest of the weight factors to a small 
value. The following tests show simulation response of the closed-loop system to a step control signal input of a c — 
5.0° over a 15 second interval. 

For the baseline comparison approximating single-objective control, the controller was derived with aircraft state 
tracking gain k qr =10, and all other k-gain scaling factors set to k — 0 .0 1 . The resulting response is shown in Figure 3. 
The baseline controller results in a CDavg=l-058e5 over the 15 second interval. The integrator winds up based on initial 
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error to the step response; since the integrator has a low weighting, the wind-up has a small effect on the response and 
results in a slow wind-down. The initial system response to the step disturbance is a strong initial control surface 
response that quickly stabilizes. 

Emphasizing drag -reduction in the control law is achieved by setting the drag weighting scale factor (kd=le4) 
higher than the remaining gains (all other gains set to 1.0), except the integrator feedback gain was set to a low value 
to be ignored by the control law (kqa=le-5). This results in the response shown Figure 4, with a drag reduction of 
93.8% from the baseline (CDavg reduced from 1.058e5 to 6.476e3). 

Including drag-reduction with state-tracking objectives (kd=kqr=le2) while minimizing other objectives (k=le-5) 
results in the response shown in Figure 5. The resulting drag CDavg=6.518e3, which is greater than an order of 
magnitude improvement over the baseline CDavg=T.058e5 (a 93.8% decrease). Figure 6 shows the result of further 
increasing the drag gain (kd=le8, kqr=le3, all other k=le-5) to the point where numerical issues begin to arise in the 
Matlab solver at the default precision. This results in a CDavg=5.531e3 (a 94.8% drag metric reduction from the 
baseline) but worse tracking performance. 

Figure 7 shows step response when integrator-performance objectives are added to the overall cost function in 
addition to drag-reduction and state-tracking (kd=kqr=kqa=le2), while minimizing other objectives (k=Te-5). Note 
that the inclusion of the integrator in this structure and in response to a large step input disturbance contributes to 
overshoot in the attack angle response. 

Emphasizing aeroelastic-control objectives, in addition to drag-reduction and state-tracking objectives, results in 
the response shown in Figure 8. The resulting drag CDavg=l.lo71e5 is slightly higher than the baseline 
CDavg=1.058e5. However, the flaps and slats have a more active response with larger required deflections. 

The results in Figure 8 show a large drag response and large actuator deflections, which can be addressed through 
increasing the drag-reduction objective gain (kd=le3) and the control cost gains (kre=krs=krf=l), while keeping 
vehicle state and aeroelastic control objectives the same (kqr=kqe=l). The results in Figure 9 show the effect of this 
modification, with a 71% relative reduction in drag (from CDavg=1.171e5 to CDavg=3.332e4) and a significantly 
reduced actuator magnitude response. 

Equally considering all objectives by setting all k gains uniformly to 1 yields the results in Figure 10. The resulting 
drag CDavg=2.013e5 has increased from the baseline due to balanced demands. An undesirable integrator wind-up 
is observed, which occurs over the first 2.5 seconds followed by a slow wind down, and the integrator response 
contributes to the overshoot and a low-amplitude low-frequency oscillatory response. 
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Figure 3. Baseline step response, state tracking objective. 
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Figure 5. Step response, emphasized drag-reduction and state-tracking objectives. 
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Figure 6. Step response, further increasing drag cost gain relative to state-tracking. 
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Figure 7. Step response with emphasized state-tracking, drag-reduction, integrator performance. 
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Figure 8. Step response with aeroelastic-control, drag, and state-tracking objectives. 
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VI. Conclusion 

The purpose of this initial study was to demonstrate feasibility of multi-objective drag-minimizing control for a 
flexible-wing aircraft with distributed actuators through a linear-quadratic control optimization framework, 
demonstrate an initial implementation based on the recently derived models, provide an assessment of expected 
performance improvements, and identify limitations and goals for future development of this system. The design and 
implementation of the control system presented in this initial study showed promising results. The multi -objective 
controller saw significant drag improvements as compared to the baseline single-objective controller. The controller 
was relatively straight-forward to derive and the solution was readily computationally-solvable. Incorporation of high- 
level control gains in the control structure allowed intuitive adjustments to be made to the closed-loop system response. 
Inclusion of the integrator structure was likewise easy to accomplish, and suggests further structural modifications to 
the controller will not be problematic. 

The results highlighted issues that will need to be addressed in further iterations of the control structure. In the 
preliminary control structure presented, response of the control surfaces are unconstrained, allowing large deflections, 
ignoring relative constraints between actuators, and allowing immediate response. The gains for tuning controller 
response provide a coarse mechanism to limit actuator response within a desired range. However, the control 
formulation should be rederived as a constrained control problem to formally address these issues and assess for 
robustness. A number of integrator issues were observed in the response of this preliminary control structure, many 
of which can be corrected through the addition of a model reference filter to produce an achievable and smoothly 
varying a c based on the plant dynamics, and further overshoot issues can be addressed through augmenting the 
feedback structure, such as adding derivative feedback paths. Once the structure of the controller has been finalized, 
an optimization methodology would facilitate gain selection for real-world implementation. The aircraft model used 
in this study is linearized; the effect of non-linear uncertainty may be of importance to consider. Direct sensing of all 
aeroelastic states will not be likely possible, and observer designs may be considered. Drag optimization in this study 
was based on an idealized drag polar, and methods for directly estimating and controlling drag may be necessary for 
achieving minimal drag on higher-fidelity models. Once these issues are addressed, future research will investigate 
integration into a larger control design framework, including extending control to handle speed and lateral modes, and 
integration with outer-loop control structures or guidance trajectories signals. Multi-objective trajectory optimization 
can be utilized to provide additional drag minimizing benefits, and trajectory-based control strategies 8 may be utilized 
to help improve tracking performance. Incorporating the wing-shape control targets found from offline optimization 
studies on the higher- fidelity numerical AASC model 5 should provide additional drag benefits. Finally, this controller 
and guidance and navigation solutions will be implemented and tested against the high-fidelity models and higher- 
fidelity vehicle simulators that are currently in development. 
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